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Abstract 

In this note four choices of outflow boundary conditions are considered for numerical conser- 
vation laws. All four methods are stable for Unear problems. For nonlinear problems examples 
are presented where either a boundary layer forms or the numerical scheme, together with the 
boundary condition is unstable due to the formation of a reflected shock. A simple heuristic 
argument is presented for determining the suitabiUty of the boundary conditions. 
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1 Introduction 


Consider the one dimensional scalar conservation law 


+ /(«)* = 0, 0 < x < 1. 


( 1 ) 


If 


a(u) = §t>0, 


then boundary data must be specified at x = 0. Since the characteristics have nonnegative slope, 
one cannot set a-priori boundary conditions at x = 1 [Lax73]. 

To compute the solution to (l) by finite differences, the interval [0, 1] is discretized into M 
pieces each of width 

1 


h = 


and a time interval At > 0 is specified. Setting 


A = 


M 


At 

Ax’ 


a numerical scheme for solving (1) is, for each m = 1,2, . . . M — 1, 

U m +I = u m — 1/2 — 

Here, 

< = «(mAz,nAf), 

the functions 


( 2 ) 


4 > m+ 1/2 — ^( W m-Jb+X> t*m-fc+2> • • • Wm+fc) 

1/2 = 

are the numerical fluxes and 

u(0 ( nAt) = g(t) 

is specified. For consistency <f> is required to satisfy 


...u) = f(u). 


Since only the role of the boundary conditions is of interest, three point schemes will be consid- 
ered. In this case, = ^(u m +i,u m ) and The Courant-Friedrichs- 

Lewy condition requires 

max |a(u)A| < 1. 
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As mentioned above, some form of artificial boundary condition is necessary at x = 1. Four 
choices of boundary conditions are discussed: 


1. Constant Extrapolation: 


u 


r*+l 

Af 


— U 


n+ 1 
Af-1> 


( 3 ) 


2. Linear Extrapolation: 


.n+l 


= 2 u 


n+l 

Af-1 


— U 


n+l 
Af — 2 > 


3. Quadratic Extrapolation: 


u 


n+l 

Af 


= 3u 


n+l 
Af — 1 


— 3u 


n+l 
Af — 2 


+ U 


n+l 

Af-3> 


( 4 ) 

( 5 ) 


4. Characteristic Extrapolation: 

1 ))- ( 6 ) 

Note that extrapolations (3), (4) and (5) are equivalent to requiring the first, second and third 
derivative, respectively, to be zero at the boundary. 

All four boundary conditions are stable for linear problems, see for example [GKS72]. Boundary 
condition (4) is of particular interest here. If ^ is a monotone function of its arguments, then 
[HHL76] showed that limit solutions of (2) are physical solutions of (1). They also showed that 
such schemes are at most first order accurate. Thus, the error using linear extrapolation is the 
same order as the error in the numerical method. In [AM81] this choice was shown to be stable 
for two dimensional explicit schemes, again for linear problems. Also, in [Pul81] this boundary 
condition was applied to the Euler equations. 

The purpose of this note is to point out that for a number of simple and well known explicit 
methods boundary conditions (4) and (5) are unsatisfactory for nonlinear problems. Examples are 
presented where, depending upon the flux function /(tt), either an artificial shock forms at the 
boundary or a reflected shock forms. 


2 Numerical Observations 


In the experiments the following initial data is used: 


ti(x,0) = u 0 (a:) = 


1 x < p, 

0 x > p. 


0 < p < 1 


( 7 ) 
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The functions /(u) are chosen so that the jump in u travels in the direction of increasing x with 
speed 

1 

8 ~ 2 ' 

Since «o > 0 the discrete solution at time t is simply 

= «o(n»Ax - 

In particular, after N time steps, where N depends on p, the numerical solution is 


u true( I > N) = 1> 0 < x < 1. 

The three choices of / are: 

/W 
/(«) 

/(«) 

Equation (9) is known as Burgers’ equation. 

Experiments were run with a number of different numerical schemes. In all of the tests, h = 1/32 
and the initial jump in u was at x = 3/4. For the numerical experiments, starting with the initial 
data in (7), 50 iterations were run and the error, 


2°’ 


= -u 


u(u--). 


( 8 ) 

( 9 ) 

( 10 ) 


€ 


n 



- «true( mAa: )) 2 


I 


1/2 


was computed. The value of A was chosen to be the smaller of the CFL condition and the stability 
bound. Note that e 50 is expected to be zero — the boundary conditions should not preclude the 
profile from passing through. 

The numerical schemes tested are: 


Lax-Friedrichs: Set 


«C +l = + -Cm) - !</-+! -/—>)• 


( 11 ) 


The Lax-Friedrichs scheme can be written in conservation form by setting 


<^( u m+ 1) u rrt) 2 ^ o \ Um y fm+1 "b ^ f m ' 


2A 
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Method 

f 

T* 

K 

«(«-*) 

Lax-Friedrichs 

0 

2.0 

oo 

MacCormick 

0 

0 

oo 

Lax-WendrofF 

0 

2.36 

oo 


Table 1: Error using Linear Extrapolation 


MacCormick’s Method: First set 


U m “ u tn ^ ( fm + 1 fm) 


and then put 




~ u m 9 (/m+ 1 fm) w(fm ~ fm-l)' 


( 12 ) 


Again (12) can be written in conservation form by setting 

1/2 2 ^ m+1 2^ m ^(/w*+l /m)) • 

Lax-Wendroff: To compute 1 by the Lax-Wendroff scheme, set 


A A 2 

U m “ U m ~ ^(/ m+1 — / m— l) + "^~[ a m+l/2 (/m+1 — /m) — a m-l/2 (/m “ /m-l)]- 


(13) 


To write the scheme in conservation form set 


^m+l/2 tyifm+l “b fm) ty ^ a m+l/2(/m+l /m)- 


In equations (11)— (13), 

/m = / ( u m) ) a m+l/2 = + 2 U m+l)> etc - 

Tables 1 and 2 display the observed value of e 50 . For constant and characteristic extrapolation 
the error was equal to zero, for each choice of /. 
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Method 

/ 

> 

2 


Lax-Friedrichs 

0 

0.0 

00 

MacCormick 

0 

00 

00 

Lax-Wendroff 

0 

00 

00 


Table 2: Error using Quadratic Extrapolation 

3 Linear and Quadratic Boundary Conditions 

From Tables 1 and 2 it is clear that when linear and quadratic extrapolation is used to compute 
the value of ujvf something fundamentally different happens when / is nonlinear. Figure 1 displays 
the solution for each of the 50 time steps for Burgers’ equation and the Lax-Friedrichs scheme. 
Note that the value of vlm-\ is fixed and that um < 0- The picture for the Lax- Wendroff scheme is 
similar. This is equivalent to overspecifying the problem data. For simplicity, only the calculations 
for linear extrapolation will be discussed. The quadratic case is similar. 

Starting with the initial data in (7), for a finite number of time steps the boundary conditions 
at x M do not affect the solution. However, when 

u Af-2 > 0 and u M-i = 0 ( 14 ) 

the boundary conditions do affect the character of the solution. The specific value of u^_ 2 will 
depend on the particular numerical fluxes used in (1). 

When updating using (2), the new value is 

u M -\ = U M -\ - M^m+1/2 - 4 > m - 1/2) 

Notice that it is the sign and magnitude of «m-i which determines how the solution will behave 
at xm- l- In successive iterations the numerical experiments indicate that the behavior of the 
solution is determined by the situation described in (14). 

The three choices of / will be considered separately. 
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3.1 / Linear 

When / is linear, in this case 

/(«) = 

the correction for all three schemes is 

A 

s \ f - 1 = 2 Um_2 > 

Here, the solution continues to travel to the right as it should. 


3.2 Burgers’ Equation 

For Burgers’ equation 



one sees that the correction, using both the Lax-Friedrichs scheme and the Lax-Wendroff scheme, 


to UM - 1 13 


SM-1 ~ 0 


since 

2) = /( U M- 2)- 

Here an artificial boundary layer forms at ®m- 1- Notice that the fact that «m-i = 0 does not 
depend on the choice of A. This same argument shows that u M -i = 0 for all successive iterations. 

As u M _ 2 ► 1 , um — ► — 1 and the total variation approaches 2. Note also that u M < 0 does not 

solve (1) since the continuous solution is positive for positive initial data. 

The situation for the MacCormick scheme is different. A direct calculation shows that the 

difference in the numerical fluxes is 

3 

«M-1 — 2 

and the solution propagates to the right. For u„_, small the speed of propagation is very slow. 
This creates a boundary layer at a = 1 for a number of time steps. Eventually, the solution does 
move off to the right. Figure 2 displays the solution for 50 time steps. 


3.3 f{u) = u{u- \) 

Finally, for 

/(u) = u(u - i) 
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f 


Boundary Condition 

\ U 


«(« - 5 ) 

«»' = 2u£i, - «K?» 

+ 

0 

- 

«s?‘ = 

+ 

+ 

+ 

1 = U M ~ M/K, +1 ) ” /( U M-l)) 

+ 

+ 

+ 


Table 3: Sign of «a/-i 


, he observations in Tables 1 and 2 indicate that the solution ^fleets off the boundary and that the 

boundary conditions are unstable. 

Calculating the correction 

«M-i = -A(^(-Ujtf_2,0) - 


for the three schemes yields: 


L-F 

A n 

(15) 

_ XJ A — 

S M-1 ~ 

M ~ 2 ' 


„Mac _ 
S M - 1 — 

- -[2 (uX,_ 2 ) 2 A 2 + (2(uJ,_ 2 ) 2 - 2(«S f _ 2 ) 8 )A + 
2 

(16) 

L-W _ 
8 M- 1 “ 

-^( 2 ( u M-2) 2a + u M-2]- 

(17) 


For all three methods the correction is less than aero and the computed solution diverges. 


4 Other Boundary Conditions 


The argument in Section 3 indicates that it is the speed and direction of propagatron of the so lutron 
It the boundary which determines what the character ofthe ariution wiU be. For the m. ml data 
L Tabie 3 summarises the sign of for the Ux-Friedrichs *ma The — for 

N/fnrCormick’s method and the Lax-Wendroff method are similar. 

MacCormick s metnoa goiution ting 

Note that it is exactly those cases where a M - 1 < 0 that interlere wi 

to the right. 
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4.1 Constant Extrapolation 

When 

U M = U M-1, 

the correction to reduces to evaluating 
At each iteration 


n 

M-2 


> ti 


n 

Af-1 


and since the functions <t> m±1/2 approximate f(u M±1/2 ), the solution will continue to 
the right. 


propagate to 


4.2 Characteristic Extrapolation 
For characteristic extrapolation, 


^ = * n M- A(/«) - /(««,_!)), (18) 

and the solution passes through the boundary as remarked in Section 2. Applying the above 

argument to compute for uj,., > 0 and uj,., = 0 one sees that equation (18) implies „» = 
0. Evaluating the correction term to update yields 

8M ~' = 


and u M _ j > 0. In this case the solution continues travelling to the right. 


5 Concluding Remarks 


In conclusion, the experimental results and the arguments in Sections 2 and 3 show that the 
oundary condition can have a dramatic effect on the solution. The stability analysis for linear 
problems does not extend to nonlinear problems and a different analysis is needed. 

To address the question of what is a ‘safe' outflow boundary condition, the numerical experi- 
ments show that erther constant or characteristic extrapolation is suitable. Unfortunately constant 
extrapolation ,s aero order accurate with respect to the meshsiae. Hence the solution at the bound- 
ary will suffer a loss of accuracy. To preserve higher order accuracy in the entire domain flctitious 

points will be needed. Characteristic extrapolation, while flrst order accurate, is more expensive 
to evaluate. 
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